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Abstract 

We estimate the two-loop perturbative corrections to zero-recoil 
matrix elements of the flavour-changing currents c^b and 07^75 b by 
calculating the terms of order nj a 2 s and substituting the dependence 
on the number of flavours by the first coefficient of the /3-function. 
Both for vector and axial vector currents, we find moderate two- 
loop corrections below 1% in magnitude. Using the Brodsky-Lepage- 
Mackenzie prescription to set the scale in the order-a s corrections in 
the MS scheme, we obtain [iy ~ 0.92- v /m;,m c and fi^ — 0.51^/m;,m c 
in the two cases. These scales are sufficiently large for perturbation 
theory to be well-behaved. The implications of our results to the 
extraction of | V&\ are briefly discussed. 
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1 Introduction 



It is now widely accepted that the measurement of the B — > D*£ v decay rate 
near zero recoil provides for the most reliable determination of the element V c b 
of the Cabibbo-Kobayashi-Maskawa matrix. The theoretical description of 
this process has a solid, model-independent foundation based on the heavy 
quark expansion, which provides a systematic expansion around the limit 
m&, m c — > oo. In this limit, QCD exhibits a spin-flavour symmetry for 
hadronic systems containing a heavy quark The symmetry-breaking 

corrections are proportional to powers of q s (toq) or 1/uiq, where we use tuq 
as a generic notation for nib and m c . These corrections can be investigated 
in a systematic way using the heavy quark effective theory (HQET) f|]-[|]. 
In particular, at zero recoil (i.e. at equal velocities of the heavy mesons) the 
representation of the flavour- changing currents in the HQET reads 

c^b^ Vv h c v Yh b v + 0(l/m 2 Q ), 

b -> VA hl-flB K + 0(l/m 2 Q ) , (1) 



where h® are the velocity-dependent heavy quark spinors of the HQET. 
Hadronic matrix elements of the effective current operators are normalized 
because of heavy quark symmetry. The coefficients r\y and t\a in ([!]) take 
into account finite renormalizations of the currents in the intermediate re- 
gion nib > > m c . They can be obtained from an on-shell matching of 
current matrix elements in QCD with the corresponding matrix elements in 
the HQET. From a measurement of the B — > D*iv decay rate near zero 
recoil, one can extract the product | V^\ T]a (1 + 5i/ m a), where 5y m 2 stands 
for non-perturbative power corrections of order (Aqcd/^q) 2 [0 0]- 

As very precise experimental data on this decay mode become available 
[§], a detailed theoretical analysis of the symmetry-breaking corrections to 
the heavy quark limit becomes increasingly important. The power corrections 
6i/ m 2 have recently been the subject of intense interest [0]-[ff2J. Here we 
shall focus on the perturbative coefficient tja and its analogue for the vector 
current, r]y. At the one-loop order, these coefficients have been known for a 



long time [13 



i 



VA =l + C F ^[<f>(z)-2], (2) 
where CV = 4/3isa colour factor, z = m c /m b , and 

6( z ) = -3 lnz-6 = + +0(ki 8 z), (3) 

vy ! 1- z 2 120 5040 V ; ' K J 

with 0(1) = 0. Using ji = ^/m h m c for the scale in the running coupling 
constant, one obtains rjv — 1.02 and tja — 0.97. Throughout this work, we 
use the input parameters m& = 4.80 GeV, m c = 1.44 GeV, and Aqcd = 0.11 
GeV in the one-loop expression for the running coupling constant in the 
MS scheme (for rif = 4). This gives a s {m b ) ~ 0.20, a s {m c ) ~ 0.29, and 
«s(v /m & m c) — 0.24. The fact that the one-loop corrections are smaller than 
the naive expectation of a s /ir ~ 10% makes one suspicious about the im- 
portance of higher-order corrections. A renormalization-group improvement 
of @ has been performed, which sums logarithms of the type (a s \nz) n , 



a s (a s \nz) n and z(a s \nz) n to all orders in perturbation theory [14]-[[TJ]. 
However, since in the case of b — > c transitions In z is not a particularly large 
parameter, one expects that the residual two-loop corrections not included 
in this procedure are as important as some of the logarithmic terms. There- 
fore, a complete two-loop calculation seems worth while. In particular, it 
would help to reduce the scale ambiguity in the above one-loop results. Un- 
fortunately, however, such a calculation appears to be rather tedious for the 
two-scale problem at hand. In this letter, we derive partial results for the 
two- loop corrections to i] V and tja, which may be used to obtain an estimate 
of the size of the full corrections. Moreover, our analysis will allow us to 
study the convergence of perturbation theory for b — > c transitions. It is thus 
of interest even beyond the two-loop order. We emphasize that our main goal 
is to investigate whether there are indications for large higher-order terms in 
the perturbative series for rj v and t/a, and not so much to obtain predictions 
for these quantities that are more precise than existing ones. To this end, it 
would be necessary to perform the complete two-loop calculations. 

Let us write the perturbative series for any one of the coefficients rjv and 
t)a in the form 

v _ 1 = f : (?M)\M i ( 4 ) 



1\ V 47T 

where a s (fi) is the running coupling constant renormalized at some scale \x 
Since rj is renormalization-group invariant, the /i-dependence on the right- 



2 



hand side must cancel between the expansion coefficients and the running 
coupling constant. It is useful to make explicit the dependence of the coeffi- 
cients t] n (fi) on the number of quark flavours. In the case at hand, the first 
dependence on rif comes at the two- loop order from diagrams containing a 
quark loop in a gluon propagator. In general, we may write 

r] n (ll) = C n , n _i(/i) Pq' 1 + C nt n- 2 (H) Po~ 2 + ■ ■ ■ + Cn,fl(n) ] 71 >1, (5) 

where Po = 11 — 1 72/ is the first coefficient of the /3-function. In particular, 
at the two-loop order we have 

v - 1 = ^ c,M + (^) 2 [cM fa + c*M] + .... (6) 

For the case of currents composed of one heavy and one light quark, it has 
been found in explicit calculations that there are large two-loop coefficients 
when the scale in the running coupling constant is chosen to be the heavy 
quark mass tjiq, and that these large coefficients are dominated by the term 
proportional to C2,i Po in @ fll9| . This empirical observation can be under- 
stood if one assumes that in these on-shell calculations the relevant scale is 
much below the "natural" scale mg, meaning that in loop diagrams virtual 
momenta below niq give a sizeable contribution. Because of the relation 

a s M = a s (fi) £ f^MY ( - ffi + . . . , (7) 

n=0 V 47T / 

using an inadequate scale can induce large higher-order coefficients c n) n-\- It 
is possible to absorb some of these large corrections by using a lower scale. 
However, in some cases this scale turns out to be too low for perturbation 
theory to be well-defined. Some of the heavy-light currents considered in 
Ref. |19| provide an example of this phenomenon. As we will discuss below, 
another example is provided by inclusive decays of hadrons containing a 



heavy quark EH]. One of our purposes here is to investigate if something 



similar happens for currents composed of two heavy quarks. 



2 Large-n/ asymptotics of perturbation the- 
ory 

The above discussion justifies that a calculation of the coefficients c„ jn _i, and 
in particular of the two-loop coefficient C2,i, is worth while. Not only can it 
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serve as an estimate of the size of the full two-loop correction, but also to 
choose an appropriate scale in the order-a s term. Technically, the coefficients 
c n)n _i can be projected out by considering the formal limit of large rif, in 
which the series (f|) takes the form 

^-1 = ^2^ [ ^ ) Cn.n-M + 0{l/n f ) 

1 oo 

= T E ln " n (^/A^cD) c n ,n-i(») + 0(l/nj) • (8) 

PO n= i 

Note that (3$ is of order n./, whereas the product j3oa s is of order n*. A 
convenient way to analyse this series is by considering its Borel transform 
with respect to ln(/i 2 /AQ CD ) [^TJ, which we define asQ 

The function -F(w, /i) can be calculated using the renormalon calculus of 



Beneke and Braun p3l E4|. In Ref. p2[, this technique has been used to 



calculate the Wilson coefficients of flavour-changing heavy quark currents 
at arbitrary velocity transfer. It is straightforward to specialize the results 
to zero recoil to obtain explicit expressions for the Borel transforms of the 
coefficients r\y and t\a- We find 



FM = C r e- (J*-Y r "g ( r 2 (1 - 2 "' ^ 

\m b m c J T(2 — u) { 2 — u 1 — . 



l-u 

z 



2(1-.),--^ + 1 ±£ 

1 + 2w 1 - z 1-^ v ; f v ' 



3C F e 



-Cu 



n \ 2u ( ^ // \ 2u ] 2A r(u) r(i - 2m) 



r(3-u) 



where again z = m c /mh. The upper (lower) sign in the first term in parenthe- 
sis refers to the vector (axial vector) current. C is a scheme- dependent con- 
stant, with C = —5/3 in the MS scheme. The scheme- and scale-dependence 
of F(u, /i) cancels against the scheme- and scale-dependence of the running 

1 This definition differs from the one adopted in Ref. (22) by a factor 1//3q. 
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coupling constant when one inverts the Borel transformation using the inte- 
gral relation 



V-l = j Q jdu (^) tt F(u )/t )+0(lft 
since the product 



A 2 
/i 2 



A 



QCD ' 



-C 



A— e 5/3 



(11) 



(12) 



is scheme- and scale- independent. 

According to @, the coefficients c n n _i can be obtained from a expansion 
of the Borel transform in powers of u. Substituting the result back into (|6|), 
we find that the //-dependence indeed cancels. At the two-loop order, we 
obtain 

2 



Vv 



VA 



1 + 



+ 



An) 



Air 



c 



[z)-2 



c 



(f)(z) Po + c 2fi (z) 



(13) 



where a s = a s (^/m^m~ c ). Since by charge conservation the vector current is 
not renormalized for z — 1, it follows that cJf (l) = 0. For 2 = m c /m& = 0.3, 
we obtain in the MS scheme (with C = —5/3) 

rjv ~ 1 + 0.236 — + (o.082 + ^c 2 %(z)) (— ) +... , 



Va 



1 



0.431 — + 

7T 



1.211 + i C ^) -f +... 



(14) 



where we use /3 = 25/3, corresponding to n/ = 4, which is appropriate for 
the intermediate region mj > \i > m c . The partial two-loop corrections that 
we have computed amount to very moderate effects, which however have 
the same sign as the one-loop corrections. Numerically, with a s ~ 0.24, 
we obtain Srjy - 5 x 10~ 4 for the corresponding contribution to rjy, and 
St)a — —7 x 10~ 3 for the contribution to tja- Thus, we find no indication 
for large two-loop corrections in the case of heavy-heavy currents. This is 
in stark contrast to the case of heavy-light currents, where the coefficient of 
the (a s /n) 2 term is typically of order 10, with ^ 02,0 of order unity [ |i~9f . 
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3 BLM scale setting 



Brodsky, Lepage and Mackenzie (BLM) have advocated to absorb vacuum 
polarization effects into the running coupling constant by choosing the scale 
in the order-a s correction so that there are no corrections of order /3q a 2 s 



in an expansion such as (y) P5| . This physically appealing scale-setting 
prescription usually results in a reasonable perturbative series. Accepting 
this point of view, one may argue that perturbation theory works well in 
cases where the BLM scale is sufficiently large, whereas it breaks down if 
this scale is too low. Recently, it has been shown that the BLM scale to be 
used in inclusive B — > X tv decays is very low, yU inc i ~ 0.07m;,, indicating 
a breakdown of perturbation theory [p0~ |. In fact, it had been noted before 



that the one-loop corrections to the inclusive decay rate exhibit a strong scale 
dependence |26| . This observation puts severe limitations on the usefulness of 
inclusive decays for the determination of | V c b\- At least, a calculation of the 
two-loop corrections is necessary before a reliable analysis can be performed. 
Fortunately, as we will now show, the situation appears to be much better 
for the exclusive decay B — > D*£ v. 

From (|6|) and fl7|), it follows that for a general perturbative series the BLM 
scale is given by 

which can be shown to be //-independent. Note that the BLM scale is not 
scheme-independent; instead, it is such that the value of (^(/^blm) is scheme- 
independent. According to ([12]), this requires that Hblm oc e°l 2 . Indeed, from 
our calculation in the previous section, we obtain 



Hv = e 3//4 e°/ 2 y/m b m, 



ha = e ^( r%j^ K /2 v^- (is) 

For z = 0.3 and in the MS scheme (with C = —5/3), this yields fiy ~ 
0.920 ^m b m c and \ia — 0.509 ^m h m c . These scales are large enough for per- 
turbation theory to be well-behaved. The corresponding scheme-independent 
values of the running coupling constant are a s (fJLv) — 0.24 and a s (fiA) — 0.30. 
Using these coupling constants instead of a s {y/m^n\.) in the one-loop ex- 
pressions (0) changes the values of r\y and t)a by 5r]y ~ 5 x 10~ 4 and 
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St)a — —9 x 10 3 . These changes are practically identical to our estimates of 
the two-loop corrections at the end of the previous section. 



4 Conclusions 

We have presented a partial calculation of the two-loop matching corrections 
to the flavour-changing currents in b — > c transitions at zero recoil. Both 
for vector and axial vector currents we find small corrections, which are 
below 1% in magnitude. Using the Brodsky-Lepage-Mackenzie prescription 
to set the scale in the order-a s corrections, we obtain a s (/iy) ~ 0.24 and 
a s (nA) — 0.30 for the relevant coupling constants in the two cases. The fact 
that the corresponding scales are sufficiently large indicates good convergence 
of perturbation theory for exclusive b —>■ c transitions. This is in contrast 
to inclusive B decays, where it was found that the appropriate scale is as 



low as 350 MeV, indicating a breakdown of perturbation theory p0 |. These 
findings support that the exclusive semileptonic decay B — > D*lv is the 
"gold-plated" mode for a precision measurement of | V c b\ [P~2 . 



It is interesting to compare our estimate of the size of the two-loop cor- 
rections with the intrinsic uncertainty in i]y and tja, which results from the 
necessity to regularize the divergent asymptotic behaviour of the perturba- 
tive series for these quantities. In fact, both r/y and tja are known to contain 
infrared renormalons, which lead to ambiguities of order (Aqcd/^q) 2 - A 



measure of the resulting intrinsic uncertainty is [22 



a 3/3 



N , 1 1 

Am ( — =f — 

m c m b 



1 2 



(17) 



where Am ~ Aqcd is the renormalon ambiguity in the pole mass of a heavy 
quark 0, |27|. As previously, the upper (lower) sign refers to the vector 
(axial vector) current. Assuming Am ~ 0.1 GeV, we obtain Ar]y ~ 0.2% 
and ArjA — 0.6%. These numbers match nicely with our estimate of the 
two-loop corrections, indicating that it is sufficient and adequate to truncate 
the perturbative series at the two-loop order. 

In conclusion, we thus believe that the existing calculations of the match- 
ing corrections for heavy quark currents, for instance the one-loop results 
r/v — 1.02 and r]A — 0.97 obtained from (fj), or the values r\y = 1.025 ±0.015 
and rjA = 0.985 ± 0.015 obtained by performing a next-to-leading order 
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renormalization-group improvement are reliable at the level of a few 
per cent. There are no indications for unusually large higher-order correc- 
tions. 
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